Search results

A finite element method is used to study the effect of flow past a circular cylinder with an integral wake splitter. A fractional step algorithm is employed to solve the Navier‐Stokes and Energy equations with a Galerkin weighted residual formulation. The vortex shedding process is simulated and the effect of splitter addition on the time period of shedding is studied at a Reynolds number of 200 and a blockage ratio of 0.25. The effect of splitter and the Strouhal number and heat transfer augmentation per unit pressure drop has been investigated.

Subcooled flow boiling phenomenon is characterized by coolant phase change in the vicinity of the heated wall. Although coolant phase change from liquid to vapour phase significantly enhances the heat transfer coefficient due to latent heat of vaporization, eventually the formed vapor bubbles may coalesce and deteriorate the heat transfer from the heated wall to the liquid phase. Due to the poor heat transfer characteristics of the vapour phase, the heat transfer rate drastically reduces when it reaches a specific value of wall heat flux. Such a threshold value is identified as critical heat flux (CHF), and the phenomenon is known as departure from nucleate boiling (DNB). An accurate prediction of CHF and its location is critical to the safe operation of nuclear reactors. Therefore, the present study aims at the prediction of DNB type CHF in a hexagonal sub-assembly.

Computational fluid dynamics (CFD) simulations are performed to predict DNB in a hexagonal sub-assembly. The methodology uses an Eulerian–Eulerian multiphase flow (EEMF) model in conjunction with multiple size group (MuSiG) model. The breakup and coalescence of vapour bubbles are accounted using a population balance approach.

Bubble departure diameter parameters in EEMF framework are recalibrated to simulate the near atmospheric pressure conditions. The predictions from the modified correlation for bubble departure diameter are found to be in good agreement against the experimental data. The simulations are further extended to investigate the influence of blockage (b) on DNB type CHF at low operating pressure conditions. Larger size vapour bubbles are observed to move away from the corner sub-channel region due to the presence of blockage. Corner sub-channels were found to be more prone to experience DNB type CHF compared to the interior and edge sub-channels.

An accurate prediction of CHF and its location is critical to the safe operation of nuclear reactors. Moreover, a wide spectrum of heat transfer equipment of engineering interest will be benefited by an accurate prediction of wall characteristics using breakup and coalescence-based models as described in the present study.

Simulations are performed to predict DNB type CHF. The EEMF and wall heat flux partition model framework coupled with the MuSiG model is novel, and a detailed variation of the coolant velocity, temperature and vapour volume fraction in a hexagonal sub-assembly was obtained. The present CFD model framework was observed to predict the onset of vapour volume fraction and DNB type CHF. Simulations are further extended to predict CHF in a hexagonal sub-assembly under the influence of blockage. For all the values of blockage, the vapour volume fraction is found to be higher in the corner region, and thus the corner sub-channel experiences CHF. Although DNB type CHF is observed in corner sub-channel, it is noticed that the presence of blockage in the interior sub-channel promotes the coolant mixing and results in higher values of CHF in the corner sub-channel.

Better membrane oxygenators need to be developed to enable efficient gas exchange between venous blood and air.

Optimal design and analysis of such devices are achieved through mathematical modeling tools such as computational fluid dynamics (CFD). In this study, a control volume-based one-dimensional (1D) sub-channel analysis code is developed to analyze the gas exchange between the hollow fiber bundle and the venous blood. DIANA computer code, which is popular with the thermal hydraulic analysis of sub-channels in nuclear reactors, was suitably modified to solve the conservation equations for the blood oxygenators. The gas exchange between the tube-side fluid and the shell-side venous blood is modeled by solving mass, momentum and species conservation equations.

Simulations using sub-channel analysis are performed for the first time. As the DIANA-based approach is well known in rod bundle heat transfer, it is applied to membrane oxygenators. After detailed validations, the artificial membrane oxygenator is analyzed for different bundle sizes (L/W) and bundle porosity (epsilon) values, and oxygen saturation levels are predicted along the bundle. The present sub-channel analysis is found to be reasonably accurate and computationally efficient when compared to conventional CFD calculations.

This approach is promising and has far-reaching ramifications to connect and extend a well-known rod bundle heat transfer algorithm to a membrane oxygenator community. As a variety of devices need to be analyzed, simplified approaches will be attractive. Although the 1D nature of the simulations facilitates handling complexity, it cannot easily compete with expensive and detailed CFD calculations.

This work has high practical value and impacts the design community directly. Detailed numerical simulations can be validated and benchmarked for future membrane oxygenator designs.

Future membrane oxygenators can be designed and analyzed easily and efficiently.

The DIANA algorithm is popularly used in sub-channel analysis codes in rod bundle heat transfer. This efficient approach is being implemented into membrane oxygenator community for the first time.

The purpose of this paper is to achieve high‐performance aerofoils that enable delayed stall conditions and achieve high lift to drag ratios.

The unsteady Reynolds averaged Navier‐Stokes equations are employed in conjunction with a shear stress transport (κ‐ω) turbulence model. A control equation is designed and implemented to determine the temporal response of the actuator. A rotating element, in the form of an actuator disc, is embedded on the leading edge of NACA 0012 aerofoil, to inject momentum into the wake region. The actuator disc is rotated at different angular speeds, for angles of attack (α) between 00 and 240.

Phenomena such as flow separation, wake vortices, delayed stall, wake control, etc. are numerically investigated by means of streamlines, streaklines, isobars, etc. Streamwise and cross‐stream forces on the aerofoil are obtained. The influence of momentum injection parameter (ξ) on the fluid flow patterns, and hence on the forces acting on the streamlined body are determined. A synchronization‐based coupling scheme is designed and implemented to achieve annihilation of wake vortices. A delayed stall angle resulted with an attendant increase in maximum lift coefficient. Due to delay and/or prevention of separation, drag coefficient is also reduced considerably, resulting in a high‐performance lifting surface.

The practicality of momentum injection principle requires both wide ranging and intensive further studies to move forward beyond the proof of concept stage.

Determination of forces and moments on an aerofoil is of vital interest in aero‐dynamic design. Perhaps, runways of the future can be shorter and/or more pay load can be carried by an aircraft, for the same stall speed.

The paper describes how a synchronization‐based coupling scheme is designed and implemented along with the RANS solver. Furthermore, it is tested to verify the dynamic adaptability of the wake vortex annihilation for NACA 0012 aerofoils.

Flow past an isolated circular cylinder and two cylinders in tandem is numerically simulated, under the influence of buoyancy aiding and opposing the flow. A modified velocity correction method is employed, which has second order accuracy in both space and time. The influence of buoyancy on the temporal fluid flow patterns is investigated, with respect to streamlines, isotherms and streaklines. Comparisons are made with respect to mean center line velocities, drag coefficients, Strouhal number and streakline patterns. Degeneration of naturally occurring Kármán vortex street into a twin eddy pattern is noticed in the Reynolds number (Re) range of 41‐200, under buoyancy aided convection. On the contrary, buoyancy opposed convection could trigger vortex shedding even at a low Re range of 20‐40, where only twin eddies are found in the natural wake. Temporal evolution of unsteady eddy patterns is visualized by means of numerical particle release (NPR). Zones of vortex shedding and twin vortices are demarcated on a plot of Richardson number against Strouhal number. Root mean square (RMS) lift coefficients (C_L,RMS) and average drag coefficient (\overline C_d) are obtained as a function of Richardson number (Ri).

The purpose of this paper is to investigate the effect of narrow gap on the fluid flow and heat transfer through an eccentric annular region is numerically. Flow through an eccentric annular geometry is a model problem of practical interest.

The approach involves standard finite volume-based SIMPLE scheme. The numerical simulations cover the practically relevant Reynolds number range of 104-106.

In the narrow gap region, temperature shoot up was observed due to flow maldistribution with an attendant reduction in the heat removal from the wall surfaces. CFD analysis is presented with the aid of, streamlines, isotherms, axial velocity contours, etc. The engineering parameters of interest such as, Nusselt number, wall shear stress, etc., is presented to study the effect of eccentricity and radius ratio.

The present investigation is a simplified model for the rod bundle heat transfer studies. However, the detailed study of sectorial mass flux distribution is a useful precursor to the thermal hydraulics of rod bundles.

For nuclear reactor fuel rods, the effect of eccentricity is going to be detrimental and might lead to the condition of critical heat flux. A thorough sub-channel analysis is very useful.

Nuclear safety standards require answers to a wide a range of what-if type hypothetical scenarios to enable preparedness. This study is a highly simplified model and a first step in that direction.

The narrow gap region has been systematically investigated for the first time. A detailed sectorial analysis reveals that, flow maldistribution and the attendant temperature shoot up in the narrow gap region is detrimental to the safe operation.

This study aims to explore an active control strategy for attenuation of in-line and transverse flow-induced vibration (FIV) of two tandem-arranged circular cylinders.

The control system is based on the rotary oscillation of cylinders around their axis, which acts according to the lift coefficient feedback signal. The fluid-solid interaction simulations are performed for two velocity ratios (V_r = 5.5 and 7.5), three spacing ratios (L/D = 3.5, 5.5 and 7.5) and three different control cases. Cases 1 and 2, respectively, deal with the effect of rotary oscillation of front and rear cylinders, while Case 3 considers the effect of applied rotary oscillation to both cylinders.

The results show that in Case 3, the FIV of both cylinders is perfectly reduced, while in Case 2, only the vibration of rear cylinder is mitigated and no change is observed in the vortex-induced vibration of front cylinder. In Case 1, by rotary oscillation of the front cylinder, depending on the reduced velocity and the spacing ratio values, the transverse oscillation amplitude of the rear cylinder suppresses, remains unchanged and even increases under certain conditions. Hence, at every spacing ratio and reduced velocity, an independent controller system for each cylinder is necessary to guarantee a perfect vibration reduction of front and rear cylinders.

The current manuscript seeks to deploy a type of active rotary oscillating (ARO) controller to attenuate the FIV of two tandem-arranged cylinders placed on elastic supports. Three different cases are considered so as to understand the interaction of these cylinders regarding the rotary oscillation.

The purpose of this study is to numerically investigate the influence of corner radius on the flow around two square cylinders in tandem arrangements at a Reynolds number of 100.

Six models of square cylinders with corner radii R/D = 0.0, 0.1, 0.2, 0.3, 0.4 and 0.5 (where R denotes the corner radius and D denotes the characteristic dimension of the body) were studied using an immersed boundary-lattice Boltzmann method, and the results were compared with those obtained using a two-dimensional unsteady finite volume method. The cylinders were mounted in a tandem configuration (1.5 ≤ L/D ≤ 10 where L denotes the in-line separation between the cylinder centers). The simulated models were quantitatively compared to the aerodynamic force coefficients and Strouhal number. Furthermore, qualitative analysis is presented in the form of flow streamlines and vorticity contours.

The R/D and L/D values were varied to observe the variation in the flow characteristics in the gap and wake regions. The numerical results revealed two different regimes over the spacing range. The drag force on the downstream cylinder was negative for all corner radii values when the cylinders were placed at L/D = 3.0 (a single-body system). Subsequently, a sudden increase was observed in the aerodynamic forces (drag and lift) when L/D increased. A different gap value was identified in the transformation from a single-body to a two-body system for different corner radii. To verify the single-body system, a simulation was carried out with a single cylinder having a longitudinal geometric dimension equal to the tandem arrangement (L/D + D). Furthermore, in a single-body regime, the total drag of a tandem cylinder was less than that of a single cylinder, thus demonstrating the benefits of using tandem structures. A significant reduction in the aerodynamic forces and drag force was achieved by rounding the sharp corners and placing the cylinders in close proximity. An appropriate configuration of the tandem cylinders with a rounded corner of R/D = 0.4 and 0.5 at L/D = 3.0 and the range is enhanced to L/D = 4.0 for 0.0 ≤ R/D < 0.4 to achieve adequate drag reduction.

To the best of the author’s knowledge, there is a paucity of studies examining the effect of corner radius on bluff bodies arranged in a tandem configuration.

The purpose of this paper is to investigate the influence of forced, in-line oscillation of a circular cylinder on an incoming incompressible flow field at different Reynolds numbers.

A space-time finite element approach is employed to model the flow around an oscillating cylinder.

The results show that two (2S), four (2P, two pair) and three vortices (P+S, one pair and one single) are shed in each cycle. In addition, a 2P _o mode is also observed, which is similar to the 2P mode but the vortices of the 2P _o mode differ in strength. The 2P mode of vortex shedding is observed along the entire wake of the flow field and 2P _o mode in the far wake. In some cases, the vortex street is transformed as it travels towards the exit to produce new patterns. One such pattern is observed for the first time in the present work, which is referred to as 2P _o * mode. The drag and lift coefficients observed are perfectly periodic at a Reynolds number of 200 and they reach a chaotic pattern as the Reynolds number is increased to a value of 350.

Originality of the paper lies in the observation of 2P vortex shedding mode or its variants in the downstream of the cylinder.

The analysis of fluid flow and thermal transport performance inside the cavity has found numerous applications in various engineering fields, such as nuclear reactors and solar collectors. Nowadays, researchers are concentrating on improving heat transfer by using ternary nanofluids. With this motivation, the present study analyzes the natural convective flow and heat transfer efficiency of ternary nanofluids in different types of porous square cavities.

The cavity inclination angle is fixed ω = 0 in case (I) and ω=π4 in case (II). The traditional fluid is water, and Fe3O4+MWCNT+Cu/H2O is treated as a working fluid. Ternary nanofluid's thermophysical properties are considered, according to the Tiwari–Das model. The marker-and-cell numerical scheme is adopted to solve the transformed dimensionless mathematical model with associated initial–boundary conditions.

The average heat transfer rate is computed for four combinations of ternary nanofluids: Fe3O4(25%)+MWCNT(25%)+Cu(50%),Fe3O4(50%)+MWCNT(25%)+Cu(25%),Fe3O4(33.3%)+MWCNT(33.3%)+Cu(33.3%) and Fe3O4(25%)+MWCNT(50%)+Cu(25%) under the influence of various physical factors such as volume fraction of nanoparticles, inclined magnetic field, cavity inclination angle, porous medium, internal heat generation/absorption and thermal radiation. The transport phenomena within the square cavity are graphically displayed via streamlines, isotherms, local and average Nusselt number profiles with adequate physical interpretations.

The purpose of this study is to determine whether the ternary nanofluids may be used to achieve the high thermal transmission in nuclear power systems, generators and electronic device applications.

The current analysis is useful to improve the thermal features of nuclear reactors, solar collectors, energy storage and hybrid fuel cells.

To the best of the authors’ knowledge, no research has been carried out related to the magneto-hydrodynamic natural convective Fe3O4+MWCNT+Cu/H2O ternary nanofluid flow and heat transmission filled in porous square cavities with an inclined cavity angle. The computational outcomes revealed that the average heat transfer depends not only on the nanoparticle’s volume concentration but also on the existence of heat source and sink.